#!/usr/bin/python3
# coding=utf-8

import numpy as np
import os

def arctan(x):
    """反正切函数：计算arctan(x)，返回弧度值（范围：-π/2到π/2）"""
    return np.arctan(x)

def gen_golden_data_arctan():
    # 配置参数（保持原有结构）
    dtype = np.float16
    input_shape = [8, 2048]  # 保持原有维度
    
    # 生成输入数据（范围设为-20到20以覆盖函数饱和区）
    # 使用log-uniform分布更合理（可选实现见注释）
    input_x = np.random.uniform(-20, 20, input_shape).astype(dtype)
    
    # 计算ArcTan（自动处理数值稳定性）
    golden = arctan(input_x).astype(dtype)
    
    # 创建目录结构（保持兼容性）
    os.makedirs("./input", exist_ok=True)
    os.makedirs("./output", exist_ok=True)
    
    # 保留tiling参数（兼容性要求）
    tiling = np.array([16384, 8], dtype=np.uint32)
    tiling.tofile("./input/input_tiling.bin")
    
    # 保存数据文件
    input_x.tofile("./input/input_x.bin")
    golden.tofile("./output/golden.bin")
    
    # 打印统计信息（调试用）
    print(f"ArcTan测试数据生成完成：")
    print(f"输入范围: [{input_x.min():.2f}, {input_x.max():.2f}]")
    print(f"输出范围: [{golden.min():.6f}, {golden.max():.6f}] (理论极值: ±{np.pi/2:.6f})")

if __name__ == "__main__":
    gen_golden_data_arctan()